Evaluation of conduction eigenchannels of an adatom probed by an STM tip 
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Ballistic conductance through a single atom adsorbed on a metallic surface and probed by a scan- 
ning tunneling microscope (STM) tip can be decomposed into eigenchannel contributions, which can 
be potentially obtained from shot noise measurements. Our density functional theory calculations 
provide evidence that transmission probabilities of these eigenchannels encode information on the 
modifications of the adatoms local density of states caused by its interaction with the STM tip. In 
the case of open shell atoms, this can be revealed in nonmonotonic behavior of the eigenchannel's 
transmissions as a function of the tip-adatom separation. 
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A deep understanding of the transport properties of 
atomic size systems is of substantial importance and has 
received considerable interest, spurred, in particular, by 
the possible applications of nanoscale conductors in fu- 
ture electronic device technologies. Thanks to advances 
in scanning tunneling microscopy, it is already possible 
not only to probe the electronic and magnetic struc- 
ture of surfaces, but also to explore conduction prop- 
erties of atomic size systems, such as one-dimensional 
wires [l[, individual organic molecules |2-6], atomic sized 
contacts [7], or even single atoms (8l-[l3|. With regard 
to the latest case, recent experimental observations by 
Berndt and coworkers revealed a qualitative differ- 

ence in the conductance behavior depending on whether 
the STM tip was approaching a clean surface or a sin- 
gle adsorbed atom. Contrary to the sudden and unpre- 
dictable jump in the conductance for the clean surface, 
a smooth and fully reversible transition from the tunnel- 
ing to the contact regime was observed for Cu(lll) and 
Ag(lll) surfaces decorated with Cu, Ag, or Co adatoms. 
Here an increased bonding of the adatom to the surface, 
due to a dipolar contribution caused by the redistribu- 
tion of the surface charge, was believed to play a crucial 
role [ll|. 

The case of magnetic adatoms deposited on nonmag- 
netic substrates is especially attractive [15[ since the 
characteristic Kondo temperature (Tk) has been ob- 
served to be sensitive to modifications of the local density 
of states (LDOS) at the adatom owing to hybridization 
of its atomic wave functions with the STM tip [9]. Here 
we offer yet another LDOS probe. The conductance of 
the atom-sized contact can be decomposed into contri- 
butions from the individual transport eigenchannels. In 
this communication we argue that the transmission prob- 
abilities of these eigenchannels, which in principle can be 
obtained from the shot noise measurements [16|, could 
serve, in addition to Tk, as a sensitive tool to probe 
the modifications in the local electronic structure of the 
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FIG. 1: A schematic STM model. The red atoms indicate 
positions where Cu/Co atoms were introduced to the system. 
The gray areas correspond to atomic planes used for the eval- 
uation of the conductance, situated four monolayers above 
the tip (Sl) and below the adsorbed atom (Sr). 



(magnetic) adatom interacting with the (spin-polarized) 
STM tip. 

We present an ab initio study of the tunneling con- 
ductance for an STM tip approaching a Cu(001) sur- 
face decorated with a single Cu or Co adatom (Fig. 1). 
Our calculations are based on density functional theory 
[17[ within the screened Korringa-Kohn-Rostoker (KKR) 
Green's function formalism [18|. The potentials were as- 
sumed to be spherically symmetric around each atom. 
Nonetheless, the full charge density, rather than its spher- 
ically symmetric part, was taken into account. The an- 
gular momentum cut-off of / max = 3 was used. Employ- 
ing an embedding technique for the Green's functions, we 
treat the STM tip and the adsorbed atom as an impurity- 
like cluster, which perturbs the otherwise perfect struc- 
ture of two semi-infinite crystalline leads separated by the 
vacuum region [19|. The zero bias conductance was cal- 
culated using the Kubo linear response theory in the for- 
mulation of Baranger and Stone [20[. Within the mixed 
site and orbital-momentum representation of the KKR 
method, pairwise contributions to the conductance have 
been summed up in real space between two atomic planes, 
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FIG. 2: Conductance as a function of the tip-adatom separa- 
tion. The inset shows a close up of the conductance for small 
tip-adatom separation including the spin-polarized contribu- 
tions. 



Sl and Sr in Fig. 1, located in the leads perpendicular 
to the z direction in which the electrical current flows 

In Fig. 2 we present results for conductance G(d) ver- 
sus tip-adatom separation d t i P - adatom computed for two 
setups: (i) a nonmagnetic case with a Cu adatom on 
Cu(001) surface approached by a Cu tip (blue points); 
and (ii) a magnetic Co atom deposited on the Cu(001) 
surface probed by a modified Cu tip with a Co apex 
atom (red points). For the latter we assumed a paral- 
lel alignment of the Co magnetic moments. As the STM 
tip is approaching the surface the tunneling conductance 
increases exponentially and thereafter smoothly evolves 
into the contact regime, similarly to the experimental 
observation jjj [l4|. However, depending on the chemi- 
cal nature of the system, with decreasing du p - adatom the 
transition from the tunneling to contact regime emerges 
in a different way. In the case of the open shell magnetic 
Co adatom (assuming we are above Tk ~ 100 K), the 
conductance experiences a strong increase in the minor- 
ity spin channel (inset of Fig. 2). From the exponential 
behavior at large distances (> 5.5 A) we have extracted 
the work functions for the systems: 5.0 eV for the Cu 
adatom, and 5.8 eV for the Co adatom, which are typical 
for metals. The first value agrees well with the tunnel- 
ing barrier height of 5.3 eV found in STM experiments 
by Berndt et al. [14]. For small tip-adatom separation 
the conductance saturates to values ~ 2.0 e 2 /h for the 
Cu system, and ~ 3.0 e 2 /h for the Co system, which are 
typical for metallic atomic point-contacts. 

According to the Landauer approach a microscopic in- 
sight into the conductance is possible by analyzing its de- 
composition into independent conduction modes (eigen- 
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FIG. 3: Transmissions of the conduction eigenchannels as a 
function of separation between the STM tip and the adatom: 
(a) a nonmagnetic Cu system; (b) majority spin channel of the 
magnetic system composed of Cu surface with a Co adatom 
probed by the Cu tip with a Co apex atom, and (c) the cor- 
responding minority spin channel. 



channels). Namely, G = GqTy[tt^] = Go^ n T n , where 
Go = 2e 2 /h is the conductance quantum and the eigen- 
channels' transmissions T n are defined as eigenvalues of 
the matrix rr^ . The matrix elements r mn are transmis- 
sion amplitudes for scattering of incoming waves m into 
outgoing waves n in the opposite lead. In the case of 
a magnetic system, the above formulation is applied to 
each spin channel separately, neglecting spin-orbit cou- 
pling. Experimental and theoretical studies 0, EH, 23] 
of atomic wires and atomic size contacts revealed a close 
relation between the conduction channels and the chemi- 
cal nature of atoms, in particular, the occupation of their 
valence atomic orbit als. Results of our calculations for 
the contact regime are in agreement with common ex- 
pectations. Namely, due to the open shell structure of 
transition metal atoms (Co in our case), one expects sev- 
eral conduction channels being partially open and con- 
ductance values above Go — a situation different from 
noble metal single atomic wires where the conductance 
is close to Go , which is usually dominated by one channel 
per spin with almost perfect transmission. 
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FIG. 4: Spin decomposed current density in the Sr plane (see Fig. 1) for majority (top row) and minority (bottom row) spin 
channels (tip-adatom separation ranging from 2.5 A to 7.0 A). White dots represent the atomic sites. 



We have performed a quantitative analysis of the eigen- 
channel contributions to conductance using a recently 
proposed approach [24], which employs a mapping 



of the transmission matrix rr^ onto a site and orbital- 
momentum KKR basis set. Within this approach we ana- 
lyze the symmetry of the solutions of the eigenvalue prob- 
lem for tt^ . We adopt a classification of eigenchannels in 
terms of the internal degrees of freedom of the selected 
atom placed at the bottleneck of the STM tunnel junc- 
tion, i.e. the adatom. Such a classification emerges natu- 
rally if one considers the wave functions of eigenchannels 
to be projected onto the adatom's valence orbitals, which 
are expected to play a decisive role in the conduction 
process [22|. With regard to the C± v symmetry of the fee 
(001) surface, one is thus able to distinguish conduction 
modes corresponding to the one-dimensional Ti (includ- 
ing s, p Zl d z 2 contributions), T3 (d x 2_ y 2) and T4 (d xy ) 
irreducible representations, and to the two-dimensional 
T 5 ({p x , d xz },{p y , d yz }) representation. 

Employing a symmetry analysis, we have found out 
that in the case of the Cu adatom probed by the Cu 
tip, Fig. 3a, only one Ti eigenchannel contributes con- 
siderably to the conductance and dominates the trans- 
mission for any dup-adatom- Other contributions to the 
total transmission are orders of magnitude smaller. Sim- 
ilarly, there exists only one dominant Ti eigenchannel 
for the majority spin states of a magnetic Co adatom 
probed by the modified spin-polarized STM tip, Fig. 3b. 
In contrast, for the Co minority spin states at smaller 
dtip- adatom almost all symmetry contributions play a 
significant role, since the partially occupied d orbitals 
start to overlap, Fig. 3c. Around dti P - a datom « 3A 
there exist five eigenchannels: two Ti channels, the dou- 
ble degenerate channel, and the T3 channel. More- 
over, and most interestingly, the transmission of these 
eigenchannels does not always change monotonically with 
dtip-adatom- In the following we show that such a behav- 
ior is ultimately related to the variations of the LDOS at 
the Co adatom induced by the interaction with the STM 
tip. 

Simultaneously to the change of the eigenchannels' 
contributions caused by the approaching tip, we also an- 



alyzed the behavior of the complementary information, 
namely a distribution of the current density in the Sr 
plane, Fig. HJ With a smaller d t ip- a datom the current 
density for the majority spin channel defocuses but oth- 
erwise, as expected, remains circularly symmetric about 
the center of Sr. This symmetry is compatible with the 
Ti representation. In contrast, the current density for the 
minority spin states changes its nature completely and 
undergoes a strong modulation reflecting the role of the d 
orbitals. The distance range for which the current density 
changes its character corresponds to d t i v - adatom ^ 3.5 A, 
where the transmission of the minority spin electrons of 
Co is dominated by the channel (Fig. [3b). 

Understanding of the above results is possible if we 
recall the projection of the eigenchannels' wave functions 
onto the adatom's orbitals. The Ti channel of Cu and of 
the Co majority states is naturally attributed to valence 
4s orbitals, while the double degenerate T 5 channel can 
be related to the valence p x and p y orbitals. Since for Cu 
and Co majority states the weight of p states at Ep is 
very low, their contribution to conductance in the contact 
regime is factor 10 smaller than Y\. 

The behavior of the Co minority spin eigenchannels as 
a function of du p - adatom (Fig. 3c) can be explained by the 
LDOS at the adatom (Fig. 5). In contrast to the major- 
ity spin d-states of Co, which are filled, the minority spin 
states are partially occupied (Fig. 5) which gives rise to a 
magnetic moment of 1.8 \±b at the Co adatom. Moreover, 
we notice that at the contact point between the adatom 
and the STM tip, d t i P - a datom — 2.5 A, each minority 
d orbital contribution is characterized by a double-peak 
structure (Fig. 5a). Because of a mutual interaction be- 
tween the Co tip and Co adatom a dimer with bonding 
and antibonding states is formed. As the STM tip is 
being retracted from the surface, the overlap between 
atomic orbitals is gradually reduced (Fig. 5b), thus the 
splitting between bonding and antibonding states dimin- 
ishes leading to the shift of minority d states towards 
Ep (Fig. 5c) and an increase in the corresponding LDOS 
(Fig. 5d-f). 

As one can see in Fig. 5a-c, the d xy states are almost 
fully occupied, and therefore their contribution to con- 
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ductance (T4 channel) is negligible. The other four mi- 
nority d orbitals are partially occupied and contribute to 
the conductance in the contact regime. The d x 2_ y 2 states 
are located in the close vicinity to Ep, giving rise to 
the T3 channel. With increasing du p - a datom the d x i_ y i 
states accumulate at Ep (Fig. 5f) and support trans- 
mission of the Ts channel until it finally starts to decay 
exponentially at larger distances. 

A similar explanation can be provided for the minority 
T5 channel of Co. In the contact regime, this channel is 
strongly influenced by the minority d xz and d yz states 
of the Co adatom. Due to the weakening of the Co-Co 
bond with increasing of the tip-adatom separation, the 
4z, d yz LDOS peak located at 0.3 eV below Ep for small 
dtip-adatom (Fig. 5a), approaches the Fermi level leading 
to an increase of the LDOS(£^f) (Fig. 5e), and there- 
fore to the increase in the conductance. The exponential 
decay later for larger d ti p- adatorn is driven by p x , p y con- 
tributions to this channel, in analogy to the case of the 
majority states. 

Finally, two minority T± channels, Ti(I) and ri(II) in 
Fig. 3c, can be interpreted as superpositions of the 4s and 
?>d z 2 states, which share the same in-plane symmetry. At 
the contact point the Ti(I) channel (filled red triangles) 
is nearly closed and thus it is projected mainly onto the 
states with d z i symmetry, whose weight at Ep is initially 
low (Fig. 5a). The ri(II) channel (open red triangles) 
has a mainly s like character and is open at the contact 
point. When the tip is retracted from the adatom, the 
antibonding d z 2 resonance approaches Ep (Fig. 5a-c,d), 
which leads to the opening of the previously closed Ti(I) 
channel. Above the crossing point, d ~ 3. OA, the two 
channels are hybrids of both orbitals, s and d z i. The 
Ti(I) channel saturates to a transmission value of almost 
1.0 due to the d z i resonance pinned to the Fermi level. 
For larger tip-adatom separations this channel acquires 
more s weight, and decays exponentially similarly to the 
majority spin T\ channel (Fig. 5b). Simultaneously, the 
ri(II) channel acquires more d z i character and decays 
much faster due to the reduced overlap between the d 
orbitals localized at the adatom and the tip. 

Summarizing, we have demonstrated that modifica- 
tions of the adatom's electronic structure due to inter- 
actions with an STM tip can result in nonmonotonic be- 
havior of eigenchannel's transmission probabilities as a 
function of the tip-adatom separation. An experimental 
confirmation of this result is expected from the analysis of 
shot noise level measurements. Moreover, the eigenchan- 
nel decomposition procedure presented here may provide 
a deeper microscopic and material-specific insight into 
the transport phenomena when applied to various STM 
related studies. 
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SPP1243, and the " Forschungsnetzwerk Nanostrukturi- 
erte Materialien" of the state Sachsen-Anhalt. 
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FIG. 5: Left: evolution of the minority spin local density of 
states (LDOS) at the Co adatom as a function of its separation 
d to the STM tip, (a) d = 2.7 A, (b) d = 3.5 A, (c) d = 4.2 A. 
The contributions of the five d orbitals are shown. Right: 
symmetry decomposed LDOS computed at the Fermi level 
(d-f), the contributions at the adatom (closed symbols) are 
compared with those at the tip apex atom (open symbols). 
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